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ABSTRACT 

In this study, the Genetic Algorithm (GA) has been used to investigate the optimum operating conditions and 
the total annual cost of heat exchangers. It is postulated to use a COM server in order to make communication easier 
between Genetic Algorithm (GA) toolbox of MATLAB and Aspen HYSYS® software. The reason of choosing the 
proposed algorithm is to improve heat exchanger performance. This has been achieved by minimizing the total annual 
cost of shell and tube heat exchanger, then implementing the output values of the decision variables into the Aspen 
EDR® software. The calculated cost of the GA - EDR case is compared with different scenario of the same shell and 
tube heat exchanger. It can be seen that the cost for the GA-EDR case is less than the cost of the studied heat exchanger 
by 37.5% and 29.2% respectively for a new and optimized cases. 
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INTRODUCTION 

Energy conservation is one of the most important research topics in the field of oil, gas, and 
petrochemical industry. Among various types of equipment used in these plants, heat exchangers have recently 
received remarkable importance. This is due to its effective conserve energy efficiency. In order to improve the 
performance of heat exchanger, various optimization techniques and softwares (Genetic algorithm and AspenTech) 
are used in this study. Accordingly Aspen HYSYS® and Aspen EDR® software’s are used for optimization and 
performance evaluation of different types of the heat exchangers. This research has focused on heat exchangers in 
the preheat train of “Basrah Refinery” to investigate the most optimum operating conditions that provide 
maximum energy saving coupled with cost reduction. Basically the preheat train is simulated and optimized using 
Aspen Hysys® software and its petroleum refining functionality, as shown in Figure 1. These optimized preheat 
exchanger then have been exported to Aspen EDR® for further rigorous simulation, sizing, and analysis. Pressure 
drops through the shell and tube sides of each heat exchanger along the preheat train have been chosen at their 
maximum values to provide more specific results. TEMA types of all heat exchangers along the preheat train have 
been used as indicated in as-built design defined available in the datasheets. Some factors from as-builts 
exchangers TEMA datasheet have been used to simulate the new heat exchanger designs, i.e. fouling resistance, 
tube length, size, and pattern configuration, allocation of hot and cold stream, and pressure drops through both 
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shell and tube sides. GAs are adaptive techniques that can be used to solve problems through investigating and optimizing 
[1]. This technique can be classified as one of the evolutionary algorithms that can rely on the imitation of the Darwin’s 
work for natural selection and survival of the strongest [2]. The basic steps of how to give a simple cycle in GA are as 
follows: 

• Creating of an initial population. 

• Estimating of each chromosome. 

• Electing the best. 

• Generating a new population of a chromosome by applying a series of operators: Selection, Crossover, and 
Mutation. 

• Selecting survivors. 

• Stopping criteria, here the cycle will be stopped when it reaches to the best solution. 

Genetic algorithm includes unique and distinguished features such as no assumptions about the function to be 
optimized are needed. Basically, GAs optimizes the tradeoff between exploring new points in the search space and 
exploiting the information discovered. Only objective function or fitness information is needed for GAs to find an optimal 
solution, this makes them competent to traditional optimization techniques that require derivatives or other auxiliary 
knowledge [3]. In the meanwhile, these features were encouraging in the application of the genetic algorithm (GA) on 
optimization and performance analysis of heat exchangers. It’s applied on one of the heat exchangers in “preheat train” 
which is called crude oil/kerosene exchanger and item number (E01). In this work, the literature review of research 
projects, which have been conducted in the field of GA, has been summarized as follows: Xie et al. (2008) applied GA to 
optimize plate-fin compact heat exchanger. Furthermore, they generalize their work and developed an optimization 
procedure to find the minimum volume and/or annual cost, and pressure drop of heat exchangers, respectively, based on 
the e-NTU and GA technique. Within their study, they presented a case study to show the optimized results by the 
proposed method that confirmed procedures 



Figure 1: Simulation of Basrah Refinery’s Preheat Train by Aspen HYSYS® 
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Applicability [4]. Rajasekara et al. (2010) applied the Modified Genetic Algorithm (MGA) to develop and test a 
model for optimization of the early design phases of shell-and-Tube Heat Exchangers. Such MGA heat exchanger models 
effectively help designers to make precise decisions at the early phases of design. This helps effectively when there is a 
shortage of available information in the design process [5]. Iqbal et al. (2014) studied a finned double-pipe heat exchanger 
with triangular fins. They aimed at finding optimal configurations of the finned annulus. They considered laminar and 
forced convection on the exchanger’s shell-side. They employed GA as the optimizer and finite element method as the 
flow solver [6]. Sadeghzadeh et al. (2015) used a multi-objective genetic algorithm to investigate the optimal design of a 
Finned Shell-and-Tube Heat Exchanger (FSTHE). Their objective functions were as: to maximize heat transfer operation, 
and minimize total cost. Throughout the study, they considered nine decision design variables, these are: tube arrangement, 
tube pitch ratio, number of tubes, tube length, tube diameter, fin height, fin thickness, number of fins per unit length and 
baffle spacing ratio [7]. Yadav et al. (2016) studied the design optimization of compact heat exchangers. They considered 
that the task of optimization, itself, may be considered as a design process in which any possible candidates will be 
analysed and evaluated based on the requirements. They concluded that GA technique may be used for optimization of 
geometrical parameters of heat exchangers to obtain optimal results [8]. 

OPTIMIZATION MODEL OF SHELL AND TUBE HEAT EXCHANGER 


The general model that was carried out on the tubular heat exchanger for optimization, has been implemented on 
the shell and tube heat exchanger. With some changes that include the following [9]: 

The same objective function for optimization is the total annual cost which is represented by the following 
equation, the utility cost for this case study will be cancelled: 


C r — A 0 KfC Ao + A 0 Ei H y Ci + A 0 E 0 H y C 0 


0 ) 


The only constraint for the objective function is derived from the overall heat balance which is related to the rate 

equation: 


q = w u c pu {t 2 - G) = w'c' p (G - t 2 ) = U 0 A 0 At lm 
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And as follows: 
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The power loss inside and outside tubes for shell and tube heat exchanger are derived as follows: 


Power Loss Inside Tubes 
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Pressure Drop Inside Tube 


A Pi = 


2fiGfLn 

PD 0 


Fanning friction factor for turbulent flow in tubes is given by: 

f _ 0.046 

" “ (DiGi/zi;) 0 - 2 


Film heat transfer coefficient inside tubes: 


(7) 
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Rearranging Equation. (9) 

(H£t 
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( 10 ) 


Also, 


Gi = 


n.DfN t 


(ID 


Substituting Equation. (10) and Equation. (7) after substituting f L , into Equation. (6) and simplifying to result: 


Ei = 


ft 3.5 B . D 1.5^. S3 (^) ' 

2 7 , 2 . 33 ^ 1.17 


1 (1.02)(0.023 ) 2 - 5 g c D 0 pfkf- 33 C^ 


( 12 ) 


Where B, is a correction factor to account for friction due to sudden contraction, sudden expansion, and reversal 
of flow direction, dimensionless. g c represents the conversion factor in Newton’s law of motion. 

The same procedure can be used to derive the power loss equation inside the tube. This power losses equation is 
then applied for the shell except that using of hydraulic equivalent diameter. The final equation is calculated based on a 
triangular tube pitch: 


V3 n 2 nD 0 


Deqv — 8 * • 


Power Loss Outside Tubes 

AP 0 G 0 S 0 


E n = 


s 0 = 


n.L.N i-D 0 p 0 

N C DcL 


(13) 

(14) 

(15) 


Where: S n = shell-side free-flow area across shell axis. 


N c = number of clearances between tubes for flow of shell-side fluid across the shell axis. 
D c = clearance between tubes to give smallest free area across shell axis. 
n b = number of baffle spaces = number of baffles + 1. 
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The Pressure Drop in the Shell Side is Calculated as 


A P n 


2BofoN r Go 


9cPo 


(16) 


Where B 0 correction factor to account for friction due to reversal of flow direction, recrossing of tubes, and 
variation in cross section, dimensionless. N r represents the number of rows of tubes across which shell fluid flows, for 
turbulent flow across tubes. 


= r) 


-0.15 


For staggered tubes. 


b 0 = 0.23 + - 


(xr-i ) 1 - 08 

Where x T = ratio of pitch normalized transverse to flow to tube diameter. 
And: 
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Where a D = 0.33 for staggered tubes. 

Also, 

_ w a .p t 
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Thus, 
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It is clear that: 

A = _5_ 
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( 21 ) 


( 22 ) 


(23) 


And: 


OPTIMIZATION STRATEGY AND PROBLEM DEFINITION 

The objective function (Ct), the total annual cost, is affected by the surface area, the power requirement for 
pumping the fluids in the inner tubes and the outer shell. In our work, the hot kerosene flows in the shell side while cold 
crude oil flows inside the tubes of the heat exchanger. Kerosene flows with the mass flowrate of 39710kg/hr at a pressure 
of 635.8kPa. The inlet and outlet temperatures of kerosene are 191 °C and 76°C, respectively. Crude oil flows with mass 
the flowrate of 397540 kg/hr at a pressure of 2016kPa. The inlet and outlet temperatures of crude oil are 25°C and 40°C, 
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respectively. At the optimum conditions, flow in the tube and shell side is turbulent. The following mentioned parameter is 
defined and used without being changed along the optimization problem [9]: 

• Installation is 15% of the purchased cost. 

• Annual fixed charges, including maintenance, are 20% of the installed cost. 

• Purchased cost of shell and tube heat exchanger with carbon steel for the outer shell and inner tubes is about 260 
$/m2. 

• Electrical cost is 0.04 $/kWh. 

• Fouling coefficient for kerosene is 2500 W/m2. K 

• Fouling coefficient for crude is 1666.67 W/m2. K 

• The exchanger operates 7000 hr/yr. 

The selected decision variables with their upper and lower bounds are as follow: 

• Number of tube passes (NP) where 1 <NP< 8 

• Tube inside diameter (Dj) where 0.005 <D,< 0.012 (m) 

• Tube outside diameter (D 0 ) where 0.008 <D 0 < 0.022 (m) 

• Number of tubes (NT) where 180 <NT< 550 

• Shell inside diameter (IDS) where 0.5 <1DS< 1.5 (m) 

• Baffle spacing (b) where 0.1 <b< 0.3 (m) 

• Number of baffles (n b ) where 1 <n b < 4 

• Tube pitch (P t ) where 0.023 <P t < 0.04 (m) 

Additional assumptions are used as follow: 

Bj = 1; B c = n b +l; F s = 1.3 and N c . N/N T =1 

Prior to use the genetic algorithm toolbox of MATLAB, genetic operators should be defined. These operators are 
defined as: 

• Number of generations = 50. 

• Population size = 50. 

• Population type = double vector. 

• Fitness scaling = Rank. 

• Selection = Tournament (Size of tournament (N) =4). 

• Crossover = Constraint dependent. 
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• Mutation = Constraint dependent. 

• Stopping criteria = Function tolerance (IE-3) 

SOLUTION METHODOLOGY FOR OPTIMIZATION OF CASE STUDY 

The optimization procedure is applied after writing a user defined function (C T ) in MATLAB as shown in the 
Figure 2. One of the questions in this work is how to calculate the physical properties for each generated solution in the 
search space. One of the used techniques is to extract the values of these properties from thermophysical database such as 
Aspen database which is attached with the aspen HYSYS® simulation software. This was achieved by establishing a 
Component Object Model (COM) server, which enabled the communication between MATLAB and Aspen HYSYS®. 
COM is a platform-independent, distributed, object-oriented system for creating binary software components that can 
interact [10]. In other words. Aspen HYSYS® and MATLAB are able to exchange data through COM. The latter transfer 
data between the two programs in binary form regardless the programming language which was used to build any one of 
the two software’s. 



Figure 2: Optimization Algorithm for Shell and Tube Heat Exchanger 
RESULTS AND DISCUSSIONS 

The shell and tube heat exchanger under study is used to heat crude oil in Crude Oil Distillation unit at the South 
Refineries Company. By using the optimization procedure, the total annual cost for shell and tube heat exchanger at new. 
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optimum and GA-EDR design, is minimized as shown in Table 1. 


Table 1: Results for Shell and Tube Heat Exchanger 
Optimization using Genetic Algorithm 


Heat Exchanger: (E01) Kerosene / Crude Oil TEMA 

Type: AES 

Item 

Design 

As-built 

New 

Optimum 

GA-EDR 

Shell-side stream: 
Kerosene 

Inlet temperature, C 

191 

191 

191 

191 

Outlet temperature, C 

76 

73.51 

73.51 

73.51 

Tube-side stream: 
Crude oil 

Inlet temperature, C 

25 

25 

25 

25 

Outlet temperature, C 

40 

39.98 

40 

39.97 

Mass flow rates, 
Kg/hr 

Shell-side stream 

39710 

39710 

39710 

39710 

Tube-side stream 

397540 

397540 

397540 

397540 

Exchanger layout 

Shell diameter, mm 

750 

750.89 

711.2 

500 

Tube Length, mm 

4500 

4500 

4267.2 

2812 

No. of tubes 

510 

485 

432 

391 

Tube diameter / pitch, mm 

20/26 

20/26 

20/26 

18/23 

Area of heat transfer, m 2 

144 

127.6 

107.6 

55.5 

Thermal analysis 

Heat duty, MM KCal/hr 

2.694 

2.626 

2.626 

3.054 

Overall heat transfer 
coefficient (service), 

Kcal/hr.m 2 . C 

217 

231.8 

274.9 

533.2 

AP SS , kg/cm 2 

0.087 

0.056 

0.038 

0.12 

AP ts , kg/cm 2 

0.884 

0.60 

0.70 

0.80 

Baffles 

Type 

Segmental 

Segmental 

Segmental 

Segmental 

Orientation 

Horizontal 

Horizontal 

Horizontal 

Horizontal 

Baffle spacing, mm 

151 

151 

203.2 

100 

% Cut 

22 

20.62 

24.15 

21 

Cost estimation 

Aspen EDR’s calculated cost, 

$ 

- 

59,969 

52,546 

36,658 


The output values for the decision variables obtained by GA optimization for shell and tube heat exchanger are 
summarizes in Table 2. 


Table 2: Optimization Results for Decision 
Variables using Genetic Algorithm 


Decision Variables 

GA Optimized, 

8 Decision Variables 

Number of tube passes x (1) 

2 

Inside diameter of tubes x (2), m 

0.012 

outside diameter of tubes x (3), m 

0.018 

Number of tubes x (4) 

391 

Inside diameter of shell x (5), m 

0.5 

Baffle spacing x (6), m 

0.1 

Number of shell baffles x (7) 

4 

Tube pitch x (8), m 

0.023 


As it can be extracted from Figure 3, the minimum total annual cost of the GA-EDR design for the shell and tube is 
achieved after 90 generations with a minimum value of 3668 USD/year. In order to test the results obtained from GA 
optimization, the values of the decision variables mentioned in Table 2 are implemented in Aspen EDR®. It should be 
noted that the cost of heat exchanger in Aspen EDR® is calculated based on the materials cost and the operating cost plus 
to that the labour cost [11]. The cost of labour can be defined as the sum of all wages paid to employees that are defined by 
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the hourly rate of work. The results from GA-EDR are compared with the as-built, new and optimized designs of the crude 
oil heater against kerosene. The results shown in Table 1 show that the cost obtained from the GA-EDR is lower than all 
the new and optimized cases by 37.5% and 29.2% respectively. 



Figure 3: Convergence Curve of GA Optimization for 
Shell and Tube Heat Exchanger 


Score Histogram 



2000 4000 6000 8000 10000 12000 14000 16000 

Score (range) 


Figure 4: Score Histogram for Shell and 
Tube Optimization Using GA 

In Figure 4, the number individuals versus their scores are illustrated showing that most of the individuals 
achieved scores around 4000 where the optimum value lies. This also confirms the steep convergence curve in Figure 3. 

CONCLUSIONS 

The optimization process can be determined in two steps, firstly, improving the design and improving operating 
conditions through providing the optimum temperature, pressure, mass flow rates, etc.. Secondly, by using the following 
optimization tools: Aspen Tech optimizer tool and MATLAB genetic algorithm toolbox. 

To optimize the performance of heat exchangers operated in the targeted industrial facilities. Aspen HYSYS® 
software has played an effective role in improving the operating conditions of the studied heat exchangers. In contrast. 
Aspen EDR® software was used to improve the mechanical and thermal design of those heat exchangers. If we have to take 
cost figures for all heat exchangers in the Crude oil Distillation unit as shown in Table 3, the economic saving 
improvement of the process heat exchangers about 40.63%. After applying the genetic algorithm, a comparison between 
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new, optimized and GA-EDR design was investigated. The cost of the GA - EDR case showed an optimization behaviour 
compared to the cost of the studied heat exchanger at a new and optimized case. 


Table 3: The Minimum Annual Cost for Shell and 
Tube Heat Exchanger 


Unit 

Total Cost, $ 

Saving Money, % 

Unoptimized 

Designs 

Optimized 

Designs 

Crude Oil Distillation 

2,470,596 

1,466,798 

40.63 
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NOMENCLATURE 
English Symbols 


A 0 

Outside surface area of heat transfer, m 2 

a„ 

Constant, a = 0.33, dimensionless 

B 

Baffle spacing, m 

b 0 

Constant in Equation. (18) for evaluating shell-side friction factor, dimensionless 


B; Correction factor to account for friction due to sudden contraction, sudden expansion, and reversal of 

flow direction, dimensionless 

B„ Correction factor to account for friction due to reversal of flow direction, recrossing of tubes, and 

variation in cross section, dimensionless 


C 

Tube clearance, m 

C P 

Heat capacity, W/kg.°C 

CA 0 

Installed cost of heat exchanger per unit of outside-tube heat transfer area, $/m 2 

Ci 

Cost of supplying 1 J to pump the fluid through the inside of the tubes, $/J 

Co 

Cost of supplying 1 J to pump the fluid through the outside of the tubes, $/J 

C T 

Total annual variable cost for heat exchanger and its operation, $/year 

D 

Diameter, m 

Ds 

Inside shell diameter, m 

I^eqv 

Hydraulic equivalent diameter based on tube pitch, m 

D c 

Clearance between tubes to give smallest free area across shell axis, m 

E 

Power loss per unit of outside-tube heat transfer area, W/m 2 

f 

Fanning friction factor, dimensionless 

F 

Correction factor, dimensionless 

F s 

Safety factor, dimensionless 

gc 

Conversion factor in Newton’s law of motion, m/s 2 

G 

Mass velocity, kg/s.m 2 

h 

Film coefficient of heat transfer, W/m 2 .K 

Hy 

Hours of operation per year, hr/year 

k 

Thermal conductivity, W/m.K 

K f 

Annual fixed charges including maintenance, expressed as a fraction of initial cost for the completely 
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installed unit, dimensionless 

L Length of one tube, m 

L t Total length of heat exchanger, m 

n Number of tube passes, dimensionless 

n b Number of baffle spaces, number of baffles + 1, dimensionless 

N c Number of clearances between tubes for flow of shell-side fluid across shell axis, dimensionless 

N r Number of rows of tubes across which shell fluid flows, dimensionless 

N t Total number of tubes in exchanger, dimensionless 

p Pressure, pa 

P, Tube pitch, m 

q Rate of heat transfer, W 

R,i Combined resistance of tube wall and scaling or dirt, m 2 .K/W 

S 0 Shell-side free-flow area across shell axis, m 2 

t Temperature of utility fluid, °C 

t' Temperature of process fluid, °C 

U Overall heat transfer coefficient, subscript d indicates that a dirt or fouling factor is included, subscript 

o indicates based on outside area and fouling factor include, W/m 2 .K 

w Mass flow, prime refers to process fluid, kg/s 

Greek Symbols 


A A t designate temperature-difference driving force, 

AP SS Shell-side Pressure drop, kg/cm 2 

AP ts Tube-side Pressure drop, kg/cm 2 

p Absolute viscosity, kg/m.s 

p Density, kg/m 3 

Yj,*F 0 Dimensional factors for evaluations of E, and E 0 

x T Ratio of pitch transverse to flow to tube diameter, dimensionless 


Subscript 

d Scaling or dirt 

i Inside pipe or tube 

lm Log mean 
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0 

Outside pipe or tube 

u 

Utility fluid 

w 

wall 


Abbreviations 


COM 

Component Object Model 

e-NTU 

Number of Transfer Units method 

EDR 

Exchanger Design and Rating 

GA 

Genetic Algorithm 

HYSYS 

Hydrocarbon Simulation system 

MGA 

Modified Genetic Algorithm 
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